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Effective magnetic penetration depth and microwave surface impedance are derived for anisotropic 
layered superconductors in the shape of spheres and long cylinders, where the external magnetic 
field is applied in the plane of the highly conducting layers to induce out-of-plane screening currents. 
The results are extended by analytic continuation to highly anisotropic conductors and to lossy 
superconductors at high frequency. The electrodynamics for the general case of a superconductor or 
metal with arbitrary anisotropy are presented. The treatment is then specialized to layered materials 
with unixaxial anisotropy, in which the penetration depth for currents fiowing perpendicular to the 
layers, Ac, is much greater than that for in- plane currents, \a. Exact solutions are found in the limit 
\a 0, and are expected to provide an accurate representation of many experiments on cuprates 
and other layered superconductors, particularly on grain-aligned powders. 
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I. INTRODUCTION 

In many superconductors of current interest the elec- 
tronic structure has strongly anisotropic character and 
the materials can be regarded as nearly two dimen- 
sional. Examples include the cuprate high tempera- 
ture superconductors,—!^ the organic superconductors,— 
and the recently discovered iron-based oxy-pnictidesi^i^ 
This reduced dimensionality gives rise to much interest- 
ing physics, such as enhanced fluctuation effectSfi and 
may even be essential to the elevated transition temper- 
atures in two dimensional materials.* However, strong 
electrical anisotropy leads to difficulties in the inter- 
pretation of electrodynamic experiments, particularly at 
RF and microwave frequencies where, in certain impor- 
tant geometries, the requirement to form closed current 
loops creates admixtures of in-plane and out-of-plane re- 
sponses. 

Finite size effects, which occur when characteristic 
electromagnetic length scales, such as penetration depth 
A and skin depth 5, become comparable to sample size, 
present an additional level of complexity in electrically 
anisostropic materials, and are often relevant in layered 
superconductors. For instance, while the in-plane pen- 
etration depth Aa in the YBa2Cu306+i, system is fairly 
short, ranging from 0.1 /xm in optimally doped material^ 
to 2 jjLm in heavily underdoped materialfi^ the out-of- 
plane penetration depth Ac ranges from 1 /xm to 100 /xm 
across the same vange^^ becoming comparable to typi- 
cal crystal sizes. Bi2Sr2CaCu208+a; is even more electri- 
cally anisotropic, with Ac ~ 100 /xm at optimal doping.^^ 
Moreover, these are low temperature values — penetra- 
tion depths grow on approach to the superconducting 
transition, crossing over to the much larger normal-state 
skin depth above Tc. 

The study of finite-size effects in isotropic supercon- 
ductors dates from the first measurements of magnetic 
penetration depth, and solutions have been given for 
a wide range of geometries J^ii^ Finite-size effects in 



anisotropic superconductors have been considered more 
recently, with various approaches taken, depending on 
the geometry of the experiment. In all cases, a com- 
ponent of the external magnetic field is applied paral- 
lel to the highly conducting layers, inducing screening 
current loops with some component flowing perpendic- 
ular to the layers, along the crystalline c axis. In the 
simplest geometries, the magnetic field is applied per- 
pendicular to the c axis, and the sample is assumed to 
be long in the field direction, eliminating demagnetiza- 
tion effects. As we will see below, this results in a two- 
dimensional screening equation that is quite tractable 
and is solved to obtain the internal magnetic field dis- 
tribution. Gough and Exon,^^ and independently Man- 
sky and co-workerSfi^ have considered such a case for an 
electrically anisotropic superconductor with rectangular 
cross-section. Their models provide realistic representa- 
tions of experiments on platelet crystals, which are usu- 
ally thin in the c direction. Porch and Waldram, working 
on grain aligned powders, were the first to discuss the 
effective penetration depth of an electrically anisotropic 
sphere jiii^ In the limit that Ac is large compared to both 
Aq and the radius of the sphere, a, screening currents are 
weak and demagnetization effects are negligible. Current 
loops, by symmetry, are then circular, and they showed 
that the internal field profile could be related to that of 
a sphere with isotropic electrodynamics and equivalent 
isotropic penetration depth A'''^ ~ "Tf^^' ^^^^ particular 
result has found much use in the analysis of experiments 
on grain-aligned powders fiiii^i^^ and is simple to apply 
because it does not depend on Aq. In the opposite limit, 
Ac ^ a. Porch and Waldram argued that the current 
loops would again be approximately circular, and that 
the problem could be mapped onto an isotropic sphere 
with A°q -Ac. 

Here we present the first complete solution to this 
problem, solving both cylindrical and spherical models 
for values of Ac ranging from zero to much greater than 
a. We confirm that the limits discussed by Waldram and 
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Porch are correct in the case of a long cyhnder. How- 
ever, we now obtain exact results for the sphere and show 
that there are significant departures from the cylindri- 
cal case, particularly in the limit Ac < a. The paper 
is organized as follows. In Sec. |TT] we derive the elec- 
tromagnetic screening equations for an anisotropic su- 
perconductor, for both 2D and 3D geometries, and then 
show how dissipation at finite frequencies can be incor- 
porated into the electrodynamics using analytic contin- 
uation. In Sec. mil we consider the problem of a long 
cylinder with anisotropic electrodynamics. This is use- 
ful both as a lead-in to the more difficult problem of the 
anisotropic sphere, and for demonstrating how experi- 
mentally accessible quantities such as surface impedance 
and penetration depth can be obtained from the inter- 
nal field distribution of the superconductor. In Sec. IIVI 
we set up and solve the problem of a spherical sample in 
which Ac ^ Xa- This is followed by a summary of our 
results in Sec. El 



II. ELECTRODYNAMICS OF ANISOTROPIC 
SUPERCONDUCTORS 



Unconventional superconductors such as the cuprates 
and organics are usually well described by London elec- 
trodynamics, as in most instances penetration depth is 
long compared to coherence length, causing the fields 
to vary slowly over the scale of a Cooper pair. In the 
London gauge, V • j4 = and the supercurrent den- 
sity js is directly related to the vector potential A by 
js = — A • A. Here A is the London parameter. In 
general A is a tensor quantity and is what character- 
izes the anisotropic electrodynamic response of layered 
superconductors. The screening equation that describes 
the Meissner state follows from combining London elec- 
trodynamics with Maxwell's equations: 

B = ^ioH = V X A (1) 
= -Vx (A-i. j,) (2) 
= -V X (A-i • V X J?) . (3) 

In the isotropic case the London parameter is a scalar, 
A, leading to the usual screening equation for magnetic 
field, 



(4) 



where A = (/ioA)^^/^ is the London penetration depth. 
In the case of anisotropic electrodynamics, and working 
in cartesian coordinates, the inverse London parameter 
is the diagonal tensor 



Ho dmg{Xl,Xl,Xl) 



Inserting this into Eq. [3] we obtain the general form for 
the anisotropic screening equation, 



{dyyHx — dxyHy) + A^ {OzzHx ~ dxzHz) 

A^ (dzzHy — dyzHz) + A^ (dxxHy — dxyHx) 

\X^ {dxxHz — dxzHx) + A^ (dyyHz — dyzHy)j 

(6) 

where, for instance, dxy denotes the partial derivative 
g"gy , and so on. In what follows we will focus on layered 
superconductors with uniaxial symmetry, in which case 
Xx ~ Xy = Xa and A^ = Ac. 

A slight modification of the approach outlined above 
allows dissipative electrodynamics at finite frequencies 
to be treated within the same framework as the static 
screening response in the Meissner state. At finite fre- 
quency, the static electromagnetic fields are replaced by 
harmonically varying, quasistatic fields, represented by 
phasor notation. For example. 



H H{t) = Re 



|jy-e''^*} 



(7) 



In this case, penetration depth A is replaced by a complex 
skin depth 5 that is determined by the complex conduc- 
tivity. 
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CT = (71 — 1(T2 = 



(8) 



In terms of 5, the surface impedance is Zg — Rs + '\Xs — 
iujfio^- In the low frequency limit, where dissipation is 
negligible, the London results are regained, with 
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-1(72 



icij/ioA^ 
iw/ioA. 



(9) 
(10) 



In the normal state, where (7s <C (7i, (5 can be directly 
related to the dc resistivity, pdc — l/f^i: 



5=il-i) 



Pdc 

2ujHo 



(11) 



(5) 



Results will be derived in the following sections in terms 
of the purely static superconductive parameters Xa and 
Ac but, in all cases, these can be generalized by analytic 
continuation using the substitution X ^ 6. The results 
for the electrically anisotropic cylinder and sphere can 
therefore be applied to a much wider range of measure- 
ments and scenarios than layered superconductors in the 
low frequency limit, including high frequency measure- 
ments on metals. In the case of the cylinder, the results 
can even be applied to magnetic materials with the sub- 
stitution /io ~* ^J'r^J■o, where fj,r is the (complex) relative 
permittivity. This extension to the magnetic case does 
not apply to the sphere, however, as the boundary match- 
ing procedure used in the solution implicitly assumes that 
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the magnetic permeabihty of the sample is equal to that 
of vacuum. Nevertheless, the method used to solve the 
spherical case could readily be extended to handle mag- 
netic materials by adapting the boundary matching pro- 
cedure. 



III. ELECTRICALLY ANISOTROPIC 
SUPERCONDUCTING CYLINDER 

Insight into the problem of the electrically anisotropic 
sphere can be obtained from solving the much simpler 
problem of an infinitely long anisotropic cylinder. In this 
case there are no demagnetizing fields, so the surface field 
is equal to the applied field Hq. Also, symmetry dictates 
that only the axial component of the magnetic field is 
nonzero. The geometry of the long rectangular platelet, 
dealt with by Gough and Exon,'^^ also has these proper- 
ties. 

We define the axis of the cylinder to lie along the a;- 
direction, and the crystal c-axis to point in the z direc- 
tion. We place the origin of the coordinate system on the 
cylinder axis. The cylinder has radius a and the location 
of its surface is given by y'^ + z'^ = a^. With this choice 
of coordinates, the superconducting layers lie parallel to 
the xy-plane. In this geometry, the screening equation, 
Eq. [HI reduces to 



8"^ 8'^ 



(12) 



We are interested in the limit of extreme anisotropy, Ac S> 
Aq, in which to good approximation we can set A^ to 
zero. In that limit the screening equation becomes one- 
dimensional: 



Hx{y,z) 



X- H 



(13) 



This limit is particularly easy to solve — magnetic flux 
penetrates the sample in planes parallel to the supercon- 
ducting layers, and is prevented from passing through the 
layers by their perfect conductivity. We can then slice the 
cylinder up into a set of independent plates, parallel to 
the superconducting layers, and solve the screening equa- 
tion separately in each one. For each plate, the process is 
identical to the problem of the infinite slab of thickness t, 
the solution of which has been tabulated by Shoenberg^H 
We substitute into the slab solution, with a slab thick- 
ness that is function of position, t{z) — 2\J d?- — z^, to 
obtain the magnetic field in the interior of the cylinder. 



Hx{y,z) 



Hq cosh — sech 
Ac 



Ac 



(14) 



Contours of contstant Hx{y,z) are plotted in Fig. [T] for 
different values of Xc/a. Despite the extreme anisotropy 
of the electrodynamics, the magnetic field distribution 
becomes isotropic in the limits Ac — > and Ac ^ a, with 




Ac = 0.3a 




Ac = a 




FIG. 1: Contours of constant field Hx {y, z) inside the cylinder, 
for difference values of Xc/a. Despite the extreme anisotropy 
of the electrodynamics, the magnetic field distribution be- 
comes isotropic in the limits Ac and Ac > a. 



the effect that screening currents will form circular loops 
in these two situations. 

One rigorous way to calculate an effective penetration 
depth, X'^^ , is to use the power-flow method of Gough and 
Exon, at finite frequency, to calculate the effective surface 
impedance, Z^^ = iw/xoA®*^, corresponding to reactive 
power flow into the sample.— For an applied field Hx (t) = 
Hq exp(iujt), Faraday's law gives the average electric field: 

<j^E-de^ iuoiio j HAy, z)dydz . (15) 

The power flowing into the cylinder per unit length is 

Pe = \Hoe-''^' j E ■ dt (16) 

- ii/oe-'"*ic^^o j H,{y,z)dydz. (17) 



The surface impedance is related to the power per unit 
1 - 

2' 



area Pa = Pi/2iia — ^Z^^H§, implying 



A 



off 



2TTaHi 



—rr I Hxiy, z)dydz . 



(18) 



We see that the effective penetration depth is related to 
the total magnetic flux inside the sample, <^Bi which can 
be partially evaluated in from the solution for H^^y, z): 



= Ho J Hxiy,z)dydz 

r 

— 4/ioi?o / Ac tanh 



(19) 
(20) 
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The effective penetration depth is then 



2'naiJLoHo 
na 



tanh 



(21) 
(22) 



This is just the average of the infinite slab expression for 
effective penetration depth, A''*^ — Atanh(i/2A)r^ over 
one quadrant, normalized to the length of arc 7ra/2. In 
the limit Ac 0, A''*^ 2X^/71 ^ 0.637Ac. We can under- 
stand this in the following way. Screening currents in this 
limit are set by the total magnetic field and are confined 
to flow parallel to the surface. We can therefore define an 
angle dependent surface impedance Zs{6) = ^J\ujiiqp{9), 
where 6 measures the angle from the y axis, with the 
superfluid response represented by an angle-dependent 
complex resistivity, 



p{d) = icj/ioA^ sin^ 9 + lujpcX^ cos^ 



(23) 



In the limit of extreme anisotropy we set Aa — 0, then 
Zs{0) — ioj/ioAcI cos0|. The effective penetration depth is 
the angle average of \{9): 



A' 



cff 



{\{e))e^\c{\coBe\)e = -\c 



(24) 



In the opposite limit, Ac S> a, the effective penetration 
depth A"'^ ^ a/2. To extract information on Ac in this 
limit it is useful to make a connection with results for the 
isotropic cylinder. The solution in the isotropic case has 
also been presented by Shoenberg,^"^ with the result that 
the internal field distribution is given by 



Mir/X) 



Jo(ia/A) 

and the effective penetration depth is given by 



A<=ff = -iA 



Ji(ia/A) 
Jo(ia/A) 



(25) 



(26) 



Here Jo{x) and Ji{x) are Bessel functions of order zero 
and one respectively. Since the internal field distribution 
of the electrically anisotropic cylinder becomes isotropic 
in the limit Ac ^ a, we can make the following argument: 
for a given value of Ac in the anisotropic cylinder, there 
is an equivalent penetration depth A'"' for the isotropic 
cylinder that results in the same internal field distribu- 
tion, and therefore the same A°*. The necessary mapping 
can be identified from series expansions of the internal 
fields. For the anisotropic cylinder. 



Hx{y,z) 



y \J o?" — Tp" 

— Hq cosh — sech 

Ac Ac 

_ TT ( 

= ^n 2A!+2A2 + 2A2 



^ " 2A| + 2Af j + 



(27) 
(28) 
(29) 



where = y^ ^ . For the isotropic cylinder, 
Jo(WA) 



Jo(ia/A) 



^°(1 + 4A^ 



1 - 



4A2 



= n 1 - 



4A2 ' 4A2 ' + - ' 



(30) 
(31) 
(32) 



where we have used Jo (a;) ~ 1 — ^2/4 + ... One expan- 
sion can be mapped on to the other with the choice 
X"^ = = O.TOTAc. To see this another way, we 

note that in this limit it is not the surface impedance 
that is angle-averaged, but the electric field. Since 
the current loops are circular, the average electric field 
{E{r)) = {p{9))i{r) is proportional to 



A2(sin2 0>+A2(cos2 0> 



2 



(33) 



The idea of an equivalent penetration depth can be 
used across the full parameter range, even when the cur- 
rent loops are no longer circular. We have shown above 
that as Ac 0, A'''^ 2Ac/7r, which is not much different 
from the finite size limit A°'^ — > Ac/\/2. In between these 
limits, we can obtain the equivalent penetration depth 
by equating the isotropic and anisotropic expressions for 
the effective penetration depth A°*^: 



.,,,Ji(ia/A<=<i) 2Ac 



-iA' 



Jo(ia/A°i) TTfl Jq 



tanh 



Ac 



-Az 



(34) 

This has been solved numerically, and the scale factor 
A'"'/Ac plotted in Fig. 3] as a function of Ac/a. 



IV. ELECTRICALLY ANISOTROPIC 
SUPERCONDUCTING SPHERE 

The general case of the electrically anisotropic sphere 
is described by the full screening equation, Eq. [SI To 
solve it, we will carry out multi-pole expansions of the 
interior and exterior magnetic fields and match them at 
the surface of the sphere. The key to making this process 
tractable is to use the somewhat unusual choice of coor- 
dinate system shown in Fig. [2l in which external field is 
applied along the x direction and the crystal c axis points 
is the z direction, with cylindrical coordinates used inside 
the sphere and spherical coordinates outside. The moti- 
vation behind this choice is that the azimuthal angle </) 
will then be the same inside and out. 



When we specialize to the limit Aa = A^; = Aj^ 
the screening equation, Eq. [51 becomes 



Hx 

Hy 



^yy y 
I dxxHy — dxyHx 



0, 



(35) 
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a- axis 
Hofx 




FIG. 2: Coordinate system for a sphere with anisotropic elec- 
trodynamics. Magnetic field Ho is applied along the x axis, 
in the plane of the superconducting layers. Spherical coor- 
dinates (r, 9, (j)) are used outside the sphere, and cylindrical 
coordinates (p, 0, z) inside. This choice of geometry greatly 
simplifies the matching of field components at the boundary, 
because the azimuthal angle is the same inside and out. 



from the radial component: 



d_ 

dp 



— ipHp) 



(39) 



The screening equation for the radial component then 
reduces to 



d^H„ idH„ 1 d^Hr,Hp_ Hp 

p^ 



tf. + I'^tlL ^ tl^ ^ tlL = tlR , (40) 

dp'^ p dp p^ dcjP' A? 



Separating variables we have Hp ~ R{p)^{ip) and 

n R" 



p 



R 



(41) 



The azimuthal component has solutions of the form 
$(0) = Acos(i^0) -I- Bsm{i'(j)), where A and B are con- 
stants and the v are positive integers. The radial equa- 
tion is 



R"{p) + -R'ip) 



P 



p- 



Rip) = . (42) 



Keeping only solutions that are regular at and even about 
the origin, the general form for the radial solution is2i 



The infinitely strong screening response implied by 
Aq — !■ means that magnetic flux cannot cross the super- 
conducting layers and H^ = inside the sphere. Com- 
bined with V • = 0, we have dyHy = —dxH^. The 
screening equation is then 



Hx 

Hy \ = X 

H, 



2 I O TT 



^yy^y 



(36) 







which we will rewrite in cylindrical coordinates. Inside 
the sphere the screening equation for the magnetic field 
decouples into separate circular slices. We use a cylin- 
drical coordinate system (p, 0, z) in which the cylinder 
axis lies along the z direction. The azimuthal angle (p 
is measured from the x-axis. The screening equation is 
V^Jf = H/\l, with Hz = 0. Expanding this out, we 
have 



(37) 



where 



1 d / dH„ 



p dp \ dp 



1 d^H„ 



(38) 



and similarly for V^H^. Again using V • = and 
Hz = 0, the azimuthal component of H can be obtained 



R{p)^J2^- 



h{p/Xc) 
P 



(43) 



where the D^, are constants and the are modified 

(or hyperbolic) Bessel functions of the first kind, of order 
V. Including the azimuthal dependence, the most general 
form for the radial component of the interior field is 



HT{p.<i^) = Y. 



Iu{p/\c) 



{Ai, cos{iy(j)) + Bi, sm{i/(j))) 



(44) 

The sm{i>(j)) terms are ruled out by the symmetry 
Hp{—y) = Hp{y). Also, as we will discuss below, only 
the v ^ I component is required to match to the exterior 
fields. In this case the interior field simplifies to 

r,int/ , /,N /i(6')a sin 6* /i(rt/Ac) cos (?!) 

W;\p, (f>,z^a cos d) = ' n/^\ ' 

p /i(asm6'/Ac) 

(45) 

where the expression has been normalized so that 
at the surface of each slice it takes the value 
H^^^{p — asin0, 0, z = acosO) — h{0). By labeling the z 
coordinate of each slice by the polar angle of its edge, we 
are anticipating that the exterior field will be expanded 
in spherical coordinates. As shown above, the azimuthal 
component of the interior field can be obtained from ra- 
dial component by differentiation to give 

H'f{p,(t>,z^ a cos 9) 

h{e)asm0 {h{p/K)+h{p/\c)) sin0 (46) 
2Ac /i(asin6'/Ac) 
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For the exterior field we use spherical coordinates 
[r,9,(j)), where the polar axis lies along the z direction 
as shown in Fig. [2l The exterior field is a solution of 
V^Jf = 0. In the space outside the sphere there is no 
current, so we can obtain the magnetic field from a scalar 
potential, ff^'^* — — V<i>m- For a field Hq applied along 
the X direction, the most general form for $m is 



$i„ = — iJo''sin6'cos(/) + ^ ^ 



CiPp {cos e)&"^'f' 



e=i m=- 



(47) 

where the are the associated Legendre polyno- 

mials. The symmetry of the problem greatly restricts 
the number of terms in the expansion. The symmetries 
of the radial component of magnetic field, Hr{x,y, z) — 
Hr{r, 9, (/)), are 



Hri-y) = +Hriy) ^ i/,.(<^) = 

cos(m(/)) but not sm{m(p) 

Hri-X) = -Hr{x) Hr{<j)) = -Hrin - 0) 

=^ m — odd 

Hr{-Z) = +Hr{z) => Hr{e) = Hr{n - 9) 

^ £ + m — even ^ £ — odd 



(48) 
(49) 
(50) 



In addition, only the m — I term is contained in the ap- 
plied field, and there is nothing in the problem to couple 
channels of different m. For this reason we can restrict 
the expansion to to = 1. The magnetic scalar potential 
is then 



rr . ^ , QP/(cOS6') C0S( 

= -Hor sm 9 cos ^2 i+i 



(51) 



odd I 



From the gradient of $ni we obtain the three components 
of the exterior magnetic field: 



dr 



Hr{r,9, 



Hq sin 9 cos cf) + 

odd i 



Ci{i+l)Pl{cos9) COS0 (52) 



^1+2 



~ r 89 

^ c A {Plicos 9)) cos cj, (53) 
Hocos9cos(f>- 22 — ^ 
odd e 



HTir,9,, 



1 



. , , QP/(cose') sin( 
^-^ sm6^r*+^ 

odd £ 



(54) 



The next step is to match the fields at the surface of the 
sphere, r = a. The azimuthal components of the interior 
and exterior fields both have a sin (j) angle dependence 



and are straightforward to match: 

H;-^r = a, 9, 0) = fff (p = a sin 9, 9, , 
qP/(cos6') 



EC£J 
si 



odd I 



sm U a 



1+2 



h{9)a sin 9 {lo{a sin 9 /Xc) + I2{asin9/ Xc)) 
2XcIi{asin9/Xc) 



(55) 
(56) 



Matching of the radial and polar components of the mag- 
netic field is more complicated. Inside the sphere, the 
infinite conductivity we have assumed for the Cu02 lay- 
ers confines the magnetic field to lie in planes parallel 
to the xy-plane, so that H^^^ = 0. However, if we also 
assume that H°^^ = Hr cos 9 — Hq sin9 = just outside 
the sphere, we run into a contradiction. Substituting the 
expansions of Hr and Hg we obtain 



Hosin9cos9+j:'-'^^^^^^^^^ 

odd e 

^ Ci^{Pli^os9))sin9 
--Hocos9sin9 - — ^ 



7^+2 



odd e 



-Ho cos 9 sin 9 (57) 
^ c£ ((£ + 1) cos 9 Pi (cos 9) - ^P/+i (cos 9)) 

odd e 



E 

odd I 



,1+2 







which has only the trivial solution q = 0, indicating that 
Hf^^ = only in the fully penetrated limit. This shows 
that for finite Ac there is an abrupt change in the di- 
rection of H on crossing the surface (see plots on the 
right-hand side of Fig. [3]), implying the presence of sur- 
face currents. This is consistent with the assumption 
of an infinite conductivity within the Cu02 layers: by 
current continuity, there will in general be in-plane com- 
ponents of the screening current that flow at the surface. 
These surface currents must have no component along 
the z axis (the crystal c axis) and therefore flow only 
in azimuthal trajectories. As such, they do not act as 
sources for H^, showing that the azimuthal component 
of the field indeed remains decoupled from Hr and Hg. 
The matching condition for Hr and Hg is the conser- 
vation of magnetic flux. This is equivalent to match- 
ing field components normal to the surface, requiring 
H^^^ {r = a,9,(t))= H'^^ {p = a sin 9,9, (j)) sin 9. The polar 
components H^^^ and do not cross the surface, so 

do not affect fiux conservation. Their difference gives the 
surface current density, jsurfacc = r x 9{Hg^*^ — H™^). 



7 










/// X \\ \ 




\\\ \ 1 




— ■ 




c-axia 11 









a z\ 



pl anp 



a 2y = plare 



FIG. 3: Lines of magnetic flux show penetration of magnetic 
field into the electrically anisotropic sphere, as a function of 
Ac/a. Figures on the left show slices through the 2 = plane, 
parallel to the highly conducting layers. Figures on the right 
show slices through the j/ = plane, in which the c axis lies. 
Note that flux lines cannot cross the layers, giving rise to 
kinks at the surface in this orientation, c-axis penetration 
depth increases from top to bottom, with Ac/a = 0, 0.2 and 1 
respectively. 



The two matching conditions at the boundary are then: 

i^osin0 + Eodd.^^^^^^^^^ = M^)Bin^, (58) 
^ qP/(cos0) 

(59) 



odd i 
2 



h{e)a sin^ e {lo{a sin 9/\c) + h {a sin 0/Ac)) 

2Ac/i(asin6'/Ac) 
h{e)saie f{x) , 



(60) 



where f{x) — x{Io{x) + l2{x)) /2Ii{x) and x = asin9/Xc 



Eliminating h{6), and defining g{x) = /(x) — 1, we obtain 

2\ CiPlicosO) 



HaP^{cos0) = 



E 

odd e 



ie + i) + 



where we have used (cos ( 
the P/(cos6') and using 



(61) 

-sinO. Projecting onto 



2i(£+l) 

Pi (cos e)Pl, (cos 0) sin 9d9 ^ , ^ ' 6e^e, 



21+1 



(62) 



we obtain 



2^(^ + 1)2 a 



odd e 



where 

{e\g-\x)\e') 



-{i'+2){£\g-\x)\£') 
(63) 



2A, 
a 



P/ (cos 61) Ji(2^) P/,(cos 



T { a sin 9 \ i j { a sin 9 \ 



2\r 



asmf 



T { a sin 9 \ 



(64) 



In practice the problem is solved numerically by truncat- 
ing the sum over t at some finite £ = ^max and solving a 
linear system for the q. The solutions converge rapidly 
with £niax, and in practice only 5 or so terms are required. 
In the limit Ac oo, g{x) ~f and all the ce = 0, cor- 
responding to full field penetration. As demagnetization 
effects are vanishingly small in this limit, the situation is 
very similar to that of the long cylinder treated above. 

Another limit is which we might expect similar be- 
haviour to the long cylinder is Ac <C a, but this turns out 
not to be the case. To examine this limit we write 

/>7r />7r 

/ P/ (cos 6')Pii (cos 6i)d6i = / P/(cos6i)(- sin6i)d6i 
Jo Jq 

L 

P/(a;)da; 

(65) 

to obtain the leading order behaviour in Xc/a: 



,+2 /<5,,i Ac 3(2^+1) 1 



The first few terms are: 



^^""12-32- 



C3 
C5 



1 97r A, 

2 ^ 32 
2l7r Ac 
1024 IT 

llTT Ac 

2048 "q" 



(66) 

(67) 

(68) 
(69) 
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FIG. 4: The scale factor between the Ac and the equivalent 
isotropic penetration depth A'^'' varies between 2/7r = 0.637 
and l/v^S — 0.707 for the electrically anisotropic cylinder and 
between 37r/16 = 0.589 and 1/^2 = 0.707 for the electrically 
anisotropic sphere. 



As Ac ^ 0, only the dipole component survives, giving 

^3^ 



= Holl- 



sm f cos ( 



H^^^ = Ha{ 1 + — I cos « cos < 



2r3 



(70) 
(71) 
(72) 



which, when transformed back into the usual spherical 
coordinate system agree with the expressions for the 
isotropic sphere in the limit A — s- 0, as tabulated by 
Shoenbergji^ However, the higher order multipole terms 
appear to linear order in Xc/a and, as we will see below, 
this leads to the nonintuitive result that the equivalent 
penetration depth of the sphere is different from that of 
the cylinder. 

The effective penetration depth is obtained by realiz- 
ing that when we carry out a penetration depth mea- 
surement, whether by microwave cavity perturbation or 
in a low frequency susceptibility experiment, we are effec- 
tively just measuring the dipole moment of the sample, 
as the higher order multipole fields fall off too rapidly 
with r to contribute to the signal. All information is 
then contained within ci(Ac), allowing us to express the 
signal as a normalized susceptibility 



X_ 

Xo 



1 



CljXc) 

97r Ac 



16 a 



3/S Ac 



(73) 
(74) 



This can be compared with the normalized susceptibility 
of the isotropic sphere, given by Shoenberg,— 



X_ 

Xo 



3A , a 

1 coth — 



3(A)2 



3— as A ^ . 
a 



(75) 
(76) 



As with the cylinder, equating the results for isotropic 
and anisotropic cases defines A"'', which will be a function 
of Xc/a. In the limit Ac — *■ 0, comparison of the leading 
terms gives A'^'^ = 37r/16 Ac = 0.589Ac. The mapping is 
plotted at other values of Ac in Fig. 21 in the form of a 
scale factor A°^/Ac. 

That the result for the sphere in the limit Ac <C a is 
different from that of the cylinder is somewhat surprising, 
as the current loops in this limit are cylindrical and flow 
in planes parallel to the y-z plane. However, as pointed 
out above, higher order harmonics are excited to linear 
order in Ac/a. This difference was also anticipated in the 
approximate treatment by Waldram and co-workers?^ 
who pointed out that the local impedance tensor on the 
surface of the sphere is a complicated function of angle, 
with principal axes that do not in general correspond 
with lines of latitude and longitude. 

Allthough we have formulated the problem with van- 
ishingly small in-plane penetration depth, a finite value 
of Xa could be introduced perturbatively, in the spirit of 
the cavity perturbation approximation used in microwave 
experimentSi ^^'^'^ As long as Xa is much less than both 
Ac and the sample radius a, the introduction of a finite 
Xa will have a negligible effect on the internal field dis- 
tribution and hence the structure of the screening cur- 
rents. Its contribution to the average surface impedance 
and effect penetration depth could, if required, be taken 
into account, for instance by carrying out a surface inte- 
gral of the angle-dependent surface impedance Zs(0, 0), 
weighted by the square of the surface magnetic field cal- 
culated in the limit Aa = 0. 

We close this section with some comments about how 
to extend the results to incorporate the high frequency 
response of metals and lossy superconductors. As ex- 
plained in Sec. |TT1 the results for the static superconduc- 
tor carry over exactly, by analytic continuation, with the 
replacement Ac — > (5c. Here Sc is the complex, c-axis skin 
depth and can be substituted directly into the expression 
for the matrix elements (^£\g~^{x)\£''), Eq.[64l The result- 
ing linear system is then solved for the complex dipole 
moment ci, which is proportional to the measured sur- 
face impedance. When considering the dissipative case, 
the scale factors plotted in Fig. |^ are no longer useful. 
Instead, the complex function ci (Sc) should be inverted 
numerically to obtain the complex c-axis conductivity di- 
rectly from the measured surface impedance. 
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V. CONCLUSIONS 

In conclusion, we have considered the problem of field 
penetration and screening in layered superconductors 
with highly anisotropic electrodynamics. Exact results 
have been derived for samples of cylindrical and spheri- 
cal shape, in geometries in which the external magnetic 
field is applied in parallel to the superconducting lay- 
ers to induce screening currents with an out-of-plane 
component. In contrast to approximations assumed in 
earlier work ,^^'^^ we find significant departures between 
the results for spheres and cylinders in the large-sample 
limit. We have argued that when the out-of-plane pene- 
tration depth is much greater than the in-plane penetra- 
tion depth, Aa can be set to zero without loss of accuracy. 
This greatly facilitates the solution of the spherical prob- 



lem, which would otherwise be intractable. A convenient 
method for conveying the results is in terms of equiva- 
lent penetration depth in an isotropic system, A'"', and 
data relating A"'^ to Ac have been tabulated for both the 
cylinder and the sphere. We have also discussed how to 
extend the results for the static superconductor to finite 
frequency measurements on lossy superconductors and 
normal metals, using analytic continuation. 
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